Unusual structural tuning of magnetism in cuprate perovskites 
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Understanding the structural underpinnings of magnetism is of great fundamental and practical 
interest. Sei_ x Te x Cu03 alloys are model systems for the study of this question, as composition- 
induced structural changes control their magnetic interactions. Our work reveals that this structural 
tuning is associated with the position of the supposedly dummy atoms Se and Te relative to the 
super-exchange (SE) Cu-O-Cu paths, and not with the SE angles as previously thought. We use 
density functional theory, tight-binding, and exact diagonalization methods to unveil the cause of 
this surprising effect and hint at new ways of engineering magnetic interactions in solids. 
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The 3d transition metal oxides are very important ma- 
terials for they have been a source of novel and intrigu- 
ing physical phenomena such as high-T c superconductiv- 
ity, colossal magneto-resistance, and magneto-electricity. 
Not surprisingly, a lot of effort is being devoted to under- 
standing the microscopic interactions that determine the 
behavior of these systems. Here we are concerned with a 
particularly important topic, namely, the structural de- 
pendence of the magnetic couplings relevant in insulators 
(e.g., direct- and super-exchange). This question already 
received a lot of attention in early studies of magnetism 
in solids 1], and a renewed interest in it is being driven 
by current intense work on magneto-electric materials. 

Sei-^Te^CuOs alloys (STCO) (3,0 are model systems 
for the study of these issues. They crystallize in a per- 
ovskite structure that is strongly distorted because both 
Te +4 and Se +4 are relatively small. Increasing x results 
in structural distortions that, in turn, switch the mag- 
netic ground state (GS) from ferromagnetic (FM) to anti- 
ferromagnetic (AFM). In Ref. Q it is proposed that the 
key structural modification is related to one of the Cu- 
O-Cu super-exchange (SE) angles present in the system, 
OL2 in Fig. [I] The change in ot2 would cause the corre- 
sponding SE coupling (J i in Fig. EJ to switch sign, thus 
transforming the GS from FM to the AFM2 spin config- 
uration of Fig. [21 This interpretation follows the spirit of 
the well-known Anderson-Goodenough-Kanemori (AGK) 
rules which discuss the SE sign and strength as a 
function of atomic species and configurations. However, 
we should note that, to the best of our knowledge, there 
is no direct experimental evidence that AFM2 is the GS 
of the Te-rich alloys, as the spin structure is yet to be 
determined by neutron scattering measurements. 

Motivated by this appealing physical picture, we de- 
cided to study the Sei-^Te^CuOa alloys using Den- 
sity Functional Theory (DFT) and complementary tight- 
binding and many-body techniques. In this Letter we re- 
port our surprising results. We find that (i) the changes 
in SE angles with x have negligible influence on the cor- 
responding couplings and (ii) what controls the magnetic 
interactions is the position of the presumedly dummy 



atoms Se and Te with respect to some Cu-O-Cu SE 
paths. Indeed, the Se/Te atoms seem to act like a valve, 
turning the AFM SE off as they approach the Cu-O-Cu 
group. Our results thus draw a picture of STCO that is 
much more subtle than that proposed in Ref. |2| . At the 
same time, they hint at new general ways of engineering 
magnetic couplings. 

The calculations were performed within the general- 
ized gradient approximation (GGA-PBE Q) to DFT. 
We primarily used the all-electron implementation in the 




FIG. 1: Left: Unit cell of ACuOs (A=Se,Te). Right top: 
symmetry-equivalent Cu-0(2)-Cu groups in ac plane. Right 
bottom: Cu-0(1)-Cu chains along b direction. Note there 
are two types of oxygens in the unit cell. Relevant structural 
parameters are defined. 



TABLE I: Structural parameters defined in Fig. Q Values 
taken from Refs. [12, [l3 • Distances in Angstroms and angles 
in degrees. Unit cell volume Q in A 3 . 



system Q a± 



OL2 



(fo/dy ^2 G?A-Q(1) 



SeCuOs 231 122.4 2.09 3.66 127.1 1.92/2.52 3.98 1.75 
TeCuOs 245 123.5 2.06 3.63 130.5 1.90/2.61 4.11 1.90 
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FIG. 2: Spin structures considered in this work. In the "FM" 
panel, only the four Cu atoms in the unit cell are shown, and 
the exchange constants are defined. 

WIEN2k package 5], with a mixed basis that includes 
augmented plane waves and local orbitals (APW+lo). 
We used the LDA+U scheme to properly treat the 3d 
electrons of Cu 0, Q • Typical cuprate values were taken 
for U (7.5 eV) and J (1.36 eV). We also used the ul- 
trasoft pseudopotential 8] implementation in the PWscf 
package , with the LDA+U approach of Ref. ^} and 
U =6 eV. The calculation conditions were converged 
to obtain exchange constants with an accuracy better 
than 1 meV. We checked that variations of 0.5-1 eV in U 
do not change our qualitative results. We double-checked 
all our results by performing both WIEN2k and PWscf 
calculations. In all cases we got full qualitative, and rea- 
sonable quantitative, agreement. 

Raw ab initio results.- We start by considering SeCuOs 
(SCO ) an d TeCu03 (TCO) in their experimental struc- 
tures [lj, LUj. Both compounds have a 20-atom unit cell 
and differ only by small variations in atomic positions 
and lattice constants. The structure is shown in Fig. [I] 
and the relevant structural data is given in Table U 

We describe the magnetic interactions by means of a 
Heisenberg Hamiltonian H = 1/2 ■ Si • Sj in which 
we include the exchange constants J\ , J 2 , and J3 defined 
in Fig. [2j ( J\ and J2 are, respectively, associated to SE 
angles ct\ and 01.2.) We compute the J's by requiring that 
this Hamiltonian reproduces, at a classical level, the en- 
ergy differences between the spin configurations in Fig. [3 
calculated from first-principles. 

Our ab initio all-electron results for SeCuOs and 
TeCuOa are given in the first two lines of Table [DJ 
In agreement with experiment, we find FM and AFM 



TABLE II: Exchange constants of Fig.|5]calculated for various 
systems (see text). Values are given in meV and the magnetic 
ground states (GS) are indicated. The results for SCO and 
TCO are consistent with high-T expansion fits of the suscep- 
tibility data in Ref. • In TCO there is competition between 
J2 and J3, which probably leads to interesting spin dynamics. 



system 


GS 


Ji 


J2 


Js 


SCO 


FM 


-4.4 


-1.3 


-0.8 


TCO 


AFM1 


6.3 


-1.5 


-0.5 


SCO/TCO-st 


AFM1 


17.7 


-2.3 


-0.6 


TCO/SCO-st 


FM 


-14.3 


1.1 


-0.7 




FIG. 3: Calculated spin-up charge densities. Panel a: Cu- 
0(2)-Cu groups in ac plane of SeCuOs (see Fig. [I] right top). 
Dashed line marks one Cu-0(2) broken bond (see text). Pan- 
els b and c: Cu-0(1)-Cu group and neighboring A cation, for 
SeCuOs and TeCuOs, respectively (see Fig.Q right bottom). 



ground states for SCO and TCO, respectively. However, 
the calculations predict that the GS of TCO is AFM1, 
and not AFM2 as proposed in Ref. |2|. Accordingly, it 
is Ji, not J2, the magnetic coupling that changes sign 
when going from SCO to TCO. In fact, even though the 
change in c\\ is around three times smaller than that in 
0L2 (see Table P), J\ varies by about 200% of its value, 
while J 2 remains almost constant. This clearly indicates 
that the SE angles have little influence in the magnetic 
couplings of these alloys. The last two lines of Table ITTI 
show the results obtained when we consider SeCuOs in 
the TeCu0 3 structure (denoted by "SCO/TCO-st") and 
viceversa. The results confirm that it is the structure, 
and not chemical differences between Se and Te, what 
determines the magnetic GS. 

We identify the causes of these results by examining 
the electronic densities that come out of the calculations. 
Figure shows the spin- up density along Cu-0(2)-Cu 
paths in the ac plane of SCO (the TCO result is essen- 
tially the same). As it is obvious, there are Cu-0(2) 
broken bonds. This is not so surprising when one notes 
that the broken-bond distance, d' 2 , is 2.52 A while the 
other Cu-0(2) distance, ^2, is only 1.92 A (see Table [[] 
the values for TCO are similar). The typical Cu-0 dis- 
tance in cuprates is 2 A, suggesting that in SCO and 
TCO the SE contribution to J 2 will be unconventional 
and weaker than usual. In fact, it is questionable that 
the above mentioned AGK rules apply in this case, and 
it seems reasonable that J 2 is largely independent of a2. 

Figures Ob and Oh show the spin-up charge density 
along the Cu-0(1)-Cu path for SCO and TCO, respec- 
tively. As far as the Cu-0 distances are concerned (see 
Table P), this SE path is similar in both systems and a 
more conventional one. However, there is a structural 
feature that makes a big difference between SCO and 
TCO, namely, the position of the neighboring A cation 
with respect to the 0(1) atom. In SCO, c?a-o(i) is 
1.75 A, while we have 1.90 A in TCO. Accordingly, as 
the density plots in Figs. 03 and Ofc suggest, the Se +4 
cation probably perturbs the 0-2p orbitals more than 
Te +4 does. One may thus hypothesize that this pertur- 
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pute the associated Wannier functions (WF's) ^ a and 
the parameters of the corresponding TB Hamiltonian: 



1.75 1.8 1.85 

Te--O(l) distance (A) 

FIG. 4: Calculated energy of spin configurations defined in 
Fig. El for TeCuOs, as a function of the Te-O(l) distance. 
FM configuration is taken as the zero of energy. Calculated 
exchange constants are given in meV. 



bation somehow disrupts the SE mechanism and renders 
a FM Ji in SCO, while regular SE results in an AFM J t 
in TCO. To check this conjecture, we calculated the mag- 
netic interactions in TCO as a function of the Te-O(l) 
distance. The results in Fig. 21 show that, indeed, when 
Te comes close enough to 0(1) (about 1.75 A), a FM GS 
results and J\ switches sign. This is very strong evidence 
that we have identified the structural feature that con- 
trols the magnetic ground state in Sei-^Te^CuOa alloys. 

Tight-binding analysis.- Great insight can be gained 
into the ultimate causes of these effects by discussing 
them in terms of the relevant electronic interactions in 
the system (i.e., hoppings, Coulomb, and exchange). 
To do so, we have implemented a simple and powerful 
scheme to compute tight-binding (TB) Hamiltonians that 
reproduce the first-principles electronic structure. In the 
following we sketch the method. 

Consider the atomic orbitals (AO's) of an isolated 
atom n. Let us bring those AO's into the crystal and 
call them </>J^, where R is a lattice vector and a stands 
for the quantum numbers {n, /,m}. The corresponding 
Bloch-like wave function is = 5^Re* kR 0]^, where k 



is in the Brillouin zone (BZ). Let ip^j and E^j denote, 
respectively, the eigenstates and eigenvalues of the Kohn- 
Sham Hamiltonian (H KS ) at k. We define 



J 



(1) 



where / is a normalization factor and the band index j 
runs over a set of bands that we are free to choose. We 
now perform the well-known Lowdin-Mattheiss transfor- 
mation to obtain orthonormal wave functions 



(2) 



j k'ol' 

where (5' k ) -1 / 2 is derived from the overlap matrix of the 
(/> k 's as described in Ref. |l^. Then, we can easily com- 



# KS ifc> 



N 



(3) 



where the C's are the bracketed coefficients in Eq. (0 
and N is the number of k-points in the BZ. 

If in the above sums we include a small number of 
bands (e.g., the nominal Cu-3<i and 0-2p bands), the re- 
sulting WF's will significantly differ from original AO's 
(j) and incorporate the effects of the surrounding lattice. 
The more additional bands we include, the more the 
WF's will resemble the isolated-atom orbitals. Note also 
that, by construction, our TB Hamiltonians can repro- 
duce the electronic band structure with arbitrary preci- 
sion, the only limitation being the spatial cutoff beyond 
which the hoppings are neglected in practice [l5|- 

We used this scheme to construct TB Hamiltonians 
for the non-spin-polarized band structures of SCO and 
TCO obtained from pseudopotential calculations. We 
considered all valence and low- lying conduction bands, as 
we were interested not only in the Cu-3d-0-2p couplings, 
but also in how those are modified by the A cations. The 
resulting WF's are very close to the isolated- atom AO's. 

The Cu-3d-0-2p hoppings along the J2 SE path are 
very asymmetric, as the two Cu-0(2) distances differ 
greatly. For the broken-bond pair, we obtain a maximum 
hopping of 0.35 eV, while we get 0.85 eV for the other 
Cu-0(2) pair. (These are results for SCO. The situation 
is similar in TCO.) This quantitatively confirms that the 
SE contribution to J2 will be relatively small. 

Regarding the other SE path, the two Cu-O(l) pairs 
are equivalent by symmetry, and the maximum C\i-3d- 
0-2p hopping is 0.77 eV in SCO and 0.87 eV in TCO. On 
the other hand, the maximum Se-4s-0(l)-2p hopping is 
4.51 eV in SCO, while for Te-5s-0(l)-2p in TCO we get 
3.65 eV. These results indicate a relatively strong Se-4s- 
0-2p interaction in SCO, which results in 3d-2p hoppings 
about 10% smaller than the corresponding ones in TCO. 
One can thus infer that the SE contribution to J\ will be 
smaller in SCO than in TCO. 

Many-body analysis.- Our first-principles TB Hamilto- 
nians can be supplemented with Coulomb and exchange 
terms to obtain realistic models of the electronic interac- 
tions. The resulting Hubbard-like Hamiltonians make it 
possible to isolate the various contributions to the mag- 
netic couplings (direct-exchange, SE), and thus identify 
the mechanisms behind them. Here we focus on the de- 
scription of the low-energy spin excitations in particular 
Cu-O-Cu groups, where the magnetic interaction is de- 
scribed by a single exchange constant J (this will be J\ or 
J2 if we consider, respectively, Cu-0(1)-Cu or Cu-0(2)- 
Cu). The many-body states are constructed by distribut- 
ing two holes among the Cu-3d and 0-2p orbitals. We di- 
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agonalize the resulting Hamiltonian matrix and compute 
J from the energy gap between the lowest-lying singlet 
and triplet states. 

We added the following terms to our first-principles TB 
Hamiltonians: on-site Coulomb interactions for both Cu- 
3d and 0-2p electrons, on-site exchange for the 0-2p elec- 
trons (which implements Hund's first rule), and inter-site 
exchange between Cu-3<i electrons (which is well-known 
to favor FM interactions). The resulting Hamiltonians 
are a simplification of that in Eq. (4) of Ref. We 
retained the above terms partly guided by numerical ev- 
idence that they have the biggest effect in the J's. We 
treated all C\i-3d and all 0-2p orbitals equally, so that 
the model has only four parameters, namely, Coulomb 
Ud and Up, and exchange K p and Kdd- 

By choosing reasonable values of these parameters, we 
were able to obtain J\ and J2 exchange constants in 
qualitative agreement with the first-principles results for 
both SCO and TCO. We employed parameters in the 
following ranges: /7 d =8-9 eV, /7 p =5-6 eV, K p =l-2 eV, 
and i^=8.5-12.5 meV. We imposed the constraint that 
Kdd is consistent with the Cu-Cu distance, i.e., we used 
larger values for smaller distances. Our analysis led us 
to two main conclusions, (i) In all the cases, the SE 
couplings per se result in AFM interactions. The inter- 
site exchange Kdd is necessary to obtain FM J's. (ii) 
The structural changes tune the magnetic couplings via 
their effect on the magnitude of the Cu-3d-0-2p hop- 
pings. Smaller hoppings result in a smaller SE singlet- 
triplet splitting and, thus, make it easier for the Kdd 
interaction to turn the coupling FM. Given the drastic 
approximations underlying our Hubbard-like Hamiltoni- 
ans, and the fact that the U and K parameters were not 
calculated ab initio, these conclusions should be taken 
with caution. Nevertheless, they are consistent with our 
body of results and make physical sense. 

In summary, we have studied the Sei-^Te^CuOs al- 
loys in which the magnetic interactions are controlled by 
composition-dependent structural changes. We employed 
three complementary approaches, i.e., LDA+U, tight- 
binding, and exact diagonalization of effective Hamilto- 
nians for finite clusters, which allowed us to study the 
magnetic interactions in great detail. We find that the 
key structural feature is the position of the supposedly 
dummy atoms Se and Te relative to Cu-O-Cu SE paths. 
As the Se/Te atom approaches the Cu-O-Cu group, it 
acts like a magnetic valve and reduces the SE contribu- 
tion to the magnetic coupling. We find SE favors anti- 
ferromagnetism and, thus, this decrease allows direct ex- 
change between Cu-3<i electrons to render a ferromag- 
netic interaction. This is quite a surprising mechanism, 
as a more conventional one related to the SE angles might 
have been expected and, in fact, was proposed in the lit- 
erature Q. On the other hand, this effect could well 
lead to novel ways of engineering magnetic couplings in 
solids. Our work shows that complex magnetic inter- 



actions may underlay seemingly simple phenomena, and 
highlights the usefulness of ab initio studies. 

It would be interesting to confirm experimentally 
the magnetic ground state of TeCu03 that we predict. 
Also, we expect that application of external pressure on 
TeCuOs could allow to tune the magnetic couplings in 
the way shown here. Finally, we are currently investigat- 
ing other perovskites (e.g., MTi03, with M=Y, La) to 
check whether similar effects occur. 

We thank Art Ramirez, Collin Broholm, and David 
Vanderbilt for useful comments. 



[1] P.W. Anderson, Solid State Phys. 14, 99 (1963). 

[2] M.A. Subramanian, A. P. Ramirez, and W.J. Marshall, 
Phys. Rev. Lett. 82, 1558 (1999). 

[3] G. Lawes, A. P. Ramirez, CM. Varma, and M.A. Subra- 
manian, Phys. Rev. Lett. 91, 257208 (2003). 

[4] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. B 
77, 3865 (1996). 

[5] See K. Schwarz, P. Blaha, and G.K.H. Madsen, Comp. 
Phys. Comm. 147, 71 (2002), and references therein. 

[6] V.I. Anisimov, J. Zaanen, O.K. Andersen, Phys. Rev. B 
44, 943 (1993). 

[7] A.I. Liechtenstein, V.I. Anisimov, J. Zaanen, Phys. Rev. 
B 52, 5467 (1995). 

[8] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[9] S. Ba roni, A. Dal Corso, S. de Gironcoli, and P. Gian- 

nozzi , |http:/ /www . pwscf . org 

[10] M. Cococcioni and S. de Gironcoli, |cond- mat/0405160 
[11] More technical details follow. Brillouin zone integrals: 
tetrahedron method, 4x4x3 grid. WIEN2k details fol- 
low. Semicore orbitals: Cu 3s and 3p, Se 3d, and Te 
4p and 4d. Muffin-tin sphere radii: i2cu=^Se=l-8 a.u. 
and Ro—IA a.u. for SeCuOs; i?cu=^Te=l-9 a.u. and 
i?o=1.6 a.u. for TeCu03. For modified structures, 
smaller values used as needed. Plane- wave cutoff K ma ^ = 
7/i^min- PWscf details follow. Valence orbitals: Te 4d, 5s 
and 5s, Se 3d, 4s and 4p, Cu 3d and 4s, and O 2s and 2p. 
Partial core corrections used. Plane-wave cutoff: 30 Ry. 
[12] K. Kohn, K. Inoue, O. Horie, and S. Akimoto, J. Sol. St. 

Chem. 18, 27 (1976). 
[13] E. Philippot and M. Maurin, Rev. Chim. Minerale 13, 
162 (1976). 

[14] L.F. Mattheiss, Phys. Rev. B 2, 3918 (1970). 

[15] Marzari and Vanderbilt [Phys. Rev. B 56, 12847 (1997)] 
used a similar scheme, with gaussians instead of AO's, 
to obtain initial guesses of their "maximally localized" 
WF's. Recently, Ku, Rosner, Pickett, and Scalettar 
[Phys. Rev. Lett. 89, 167204 (2002)] employed such WF's 
to derive first-principles many-body Hamiltonians. Let us 
add that we obtain fairly localized WF's when a small 
number of bands (e.g., the nominal Cu-3d and 0-2p 
bands) are used in the calculation (Eq. O). Moreover, 
this construction method allows to extract the bands 
of interest in entangled-band situations and control the 
transformation properties of the resulting WF's. 

[16] T. Yildirim, A.B. Harris, A. Aharony, and O. Entin- 
Wohlman, Phys. Rev. B 52, 10239 (1995). 



